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Abstract 

The two-parameter distribution known as exponential-Poisson (EP) dis- 
tribution, which has decreasing failure rate, was introduced by Kus (2007). 
In this paper we generahze the EP distribution and show that the failure rate 
of the new distribution can be decreasing or increasing. The failure rate 
can also be upside-down bathtub shaped. A comprehensive mathematical 
treatment of the new distribution is provided. We provide closed-form ex- 
pressions for the density, cumulative distribution, survival and failure rate 
functions; we also obtain the density of the ith order statistic. We derive the 
rth raw moment of the new distribution and also the moments of order statis- 
tics. Moreover, we discuss estimation by maximum likelihood and obtain an 
expression for Fisher's information matrix. Furthermore, expressions for the 
Renyi and Shannon entropies are given and estimation of the stress-strength 
parameter is discussed. Applications using two real data sets are presented. 

keywords: Exponential-Poisson distribution, Generalized exponential-Poisson 
distribution, Failure rate, Order statistics, Fisher's information matrix. 
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1 Introduction 



Adamidis and Loukas (1998) introduced a distribution with decreasing failure 
rate. This distribution is known as exponential-geometric distribution and is ob- 
tained by compouding an exponential distribution with a geometric distribution. 
In the same fashion, Kus (2007) introduced a two-parameter distribution known as 
exponential-Poisson (EP) distribution, which has decreasing failure rate, by com- 
pouding an exponential distribution with a Poisson distribution. Proschan (1963) 
proved that the decreasing failure rate property is inherent to mixtures of distribu- 
tions with constant failure rates. The cumulative distribution function (cdf) of the 
EP distribution is given by 



x,/3,A> 0. 

Gupta and Kundu (1999) introduced the generalized exponential (GE) distri- 
bution (also known as exponentiated exponential distribution) whose the cdf is 



Following their approach, in this paper we generalize the EP distribution by 
considering the following cdf: 



a > 0. A random variable X with cdf ([T]) is said to follow a generalized exponential- 
Poisson (GEP) distribution (or exponentiated exponential-Poisson distribution). 

Other distributions obtained using this approach have been introduced and 
studied in the literature. The exponentiated WeibuU (EW) distribution introduced 
by Mudholkar et al. (1993) generalizes the GE distribution and was also consid- 
ered by Mudholkar et al. (1995), Mudholkar and Hutson (1996) and Nassar and 
Eissa (2003). Generalizations of the gamma, WeibuU, Gumbel and Frechet distri- 
butions are given by exponentiated type distributions: the exponentiated gamma, 
exponentiated WeibuU, exponentiated Gumbel and exponentiated Frechet distri- 
butions, respectively; see Nadarajah and Kotz (2006). 

Generalizations of exponentiated type distributions can be obtained from the 
class of generalized beta distributions. This class has been receiving increased 
attention recently, in particular after the works of Eugene et al. (2002) and Jones 



,-A+A exp(-/ix) 



Gep(x) = 



GGE(x) = (l-e-''T, x,A,a>0. 




(1) 
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(2004) . Eugene et al. (2002) introduced what is known as the beta normal (BN) 
distribution; the exponentiated normal distribution is a particular case of the BN 
distribution. Nararajah and Kotz (2004) introduced the beta Gumbel distribution 
which generalizes the exponentiated Gumbel distribution. Nadarajah and Kotz 

(2005) also introduced the beta exponential (BE) distribution which has the GE 
distribution as a particular case. A generalization of the exponentiated Frechet 
distribution was obtained by Nadarajah and Gupta (2004): the beta Frechet distri- 
bution. 

The remainder of our paper is organized as follows. In Section [2l we give the 
density and failure rate functions of the GEP distribution. Furthermore, we show 
that the failure rate can be decreasing or decreasing; it can also be upside-down 
bathtub. Additionally, we derive the probability density function of the ith order 
statistic. Expressions for the rth raw moments of the GEP distribution and of 
the rth order statistic are given in Section [31 We discuss estimation by maximum 
likelihood and provide an expression for Fisher's information matrix in SectionlH 
In Section|5]we give expressions for the Renyi and Shannon entropies. We discuss 
estimation of the stress-strength parameter in Section |6l Applications are given in 
the Section|71 Finally, Section[F]concludes the paper. 



2 Density, failure rate functions and order statistic 

Let X be a random variable with GEP distribution with parameters A, /3 and a, 
i.e., X ~ GEP(A,/3, a). The cdf of X is given in ([T]) and the associated probability 
density function (pdf) is 

f(x) = fl - g-'l+'ie''P(-y8-«-)|0'-lg-'l-y8-v+'le''P(--e^) (T) 

(1 ~ e~'*)" ' V / 

X, A,/3, a > 0. 

The EP distribution is a particular case of the GEP distribution corresponding 
a = I. When a = I and taking A ^ 0^,it can be shown that X ^ Y, where 
denotes convergence in distribution and 7 is a random variable with exponential 
distribution with parameter /3. 

For |z| < 1 and y > 0, we have that 

2 Izim. ,3, 
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If y is integer, the sum in ([3]) stops at y - 1 . Let a > be non-integer. Plugging 
© into ©, we obtain 



fix) 



Z 



r(a-j)(j + l)! 



(4) 



where /£p(-, + is the pdf of a random variable with EP distribution with 
parameters X{j +1) and jS. When a > is integer, the index j in the summation 
in dU ranges up to a - 1 . We then have that the the GEP pdf can be written as 
a linear combination of EP pdfs. Figure [U shows different GEP densities, i.e., it 
contains plots of 

X=\ and p=2 



a=0.1 
0=1.0 
a=5.0 
a=10 
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Figure 1: GEP densities. 
The survival function and the failure rate of the GEP distribution are given by 



s{x) = 1 



1 _ ^-A+Aexp{-/ix)Y 



, x>0, 



and 



h(x) = aA/3- 



1 - e-'^ 

^ _ ^-A+Aexp(-fSx)'^a-\ ^-A-/Sx+Aexp(-/3x) 



, X > 0, 



(5) 



respectively. 

We shall now show that the failure rate of the GEP distribution can be de- 
creasing or increasing depending on the parameter values. Define the function 
ri{x) = -f'{x)/f{x), where /' denotes the first derivative of /. It is straightfoward 
to show that 

r]{x) = -Aj3(a - \)e~\\ - e-'<+^exp(-/^^-)j-ig-/i.v+iexp(-/^x) ^^^.^^-Z^x ^ 

and 

ri'(x) = AYia - l)e-2^{l - exp(-Av)}-2g-2/3^+2iexp(-/3^) ^ ^^2(^ _ ^^^~A 

X {1 - g-''+-1exp(-/3A:)|-l| _|_ ^^-/3j:|g-/3.Y+iexp(-/3jc) _ ^'^q-^^ 

If < a < 1 , then T]'{x) < Vjc> 0. It follows from Theorem (b) of Glaser ( 1 980) 
that the failure rate is decreasing. 
Using the inequalities 

Q-A+A exp(-y3.Y) 2 

> - — - and 1 < 1 + Ae-l^' 



\ _ g-/l+iexp(-jSx) gA _ Y 

we obtain 

Therefore, if a > e'^, then t]'{x) > Vx > 0. Hence, using Theorem (b) of Glaser 
(1980), we conclude that the failure rate is increasing. Figure [21 shows plots of 
the failure rate given in ([5]). We note that, as mentioned earlier, these rates can be 
increasing or decreasing. Furthermore, we see that failure rate can be upside-down 
bathtub. 

The pdf of the /th order statistic for a random sample X\,. . .,Xn from the 
GEP distribution is given by 



and then 



fi:n{x) = -— — -fix)Fixy~'{l - Fix)}"-', 

(i- l)!(n-z)! 



f.^^^jc) = ^^^{1 - e-^+'^-P(-/^'^)}--i 



5 



;.=0.5 and |i=1 



?.=5 and p=1 




n \ \ I I I ^ I \ \ n 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 0.0 0.5 1.0 1.5 2.0 



Figure 2: GEP failure rates. 

where /(■) and F(-) are pdf and cdf of the GEP distribution, respectively. 
An useful alternative expression for the pdf of the zth order statistc is 

fi:n(x) = ——- ttV^ , ]fGEp(x;A,fi,a(i + k)), (6) 

where /gepC', A,/3, o;{i+ j)) is the pdf of a random variable having GEP distribution 
with parameters A, /3 and a{i + j). The pdf of the zth order statistic for a random 
sample of the EP distribution is obtained by setting a = 1 in Q. Expressions dH) 
and © play an important role in the derivation of the main properties of the GEP 
distribution. 

3 Quantiles and moments 

In this section we shall present closed-form expressions for the GEP quantiles and 
also for the raw moments of the distribution order statistics. The ^th quantile Xq 
of the GEP distribution can be obtained from ([T]) as 

X, = -/3-' log[l + log{l - q''"(l - e-')}] . 

In particular, the distribution median is 

xo.5 = log [l + A-' log{l - 2-"%l - e--*)}] . 
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We shall now provide a closed-form expression for the rth raw moment of the 
GEP distribution. Kus (2007) showed that the rth moment of a random variable Y 
with EP distribution with parameters A and /3 can be expressed as 

E^^') = 7T^^^^'-+i.^+i([l' • . . , 1], [2, . . . , 2], A), (7) 

where 

^^ A'nUnni + k)Y(nd-' 

with n = [ni, . . .,np] and d = [di, . . . ,dq]. The latter is the known Barnes ex- 
tended hypergeometric function, which can be easily evaluated using, e.g., the 
computer algebra systems Maple and Mathematica. 

Let X ~ GEP{A,/3, a). If a > is non-integer, using (Hj) and (|7]) we obtain that 
the rth raw moment of X can be written as 

= n Ajnr E -FT TT^r.u.iiU, . . . , 1], [2, . . . , 2], + D); (8) 

(1 -e '^r^'^ r(a-j)jl 

when a > is integer the above summation stops at a - 1. 

Let Xi-„ be the rth order statistic of a random sample from the GEP distribution. 
From we have that the rth raw moment of Xi.„ for a > non-integer can be 
written as 

nlrlA ^ P^ (-iy+*^r(a(/ + k) + l)e-'^'j+^^ 



E(x,r„) = — — y y 



X I" ^ 'ji^.+i,.+i([l, . . . , 1], [2, . . . , 2], AU + D); (9) 

when a > is integer the index j in the above summation only ranges up tp 
a - I. It is noteworthy that one can obtain the moments of the order statistics of 
the EP distribution by setting c = 1 in ©, a result which was not, to the best of 
our knowledge, available in the literature. That is, the above result yields as a by 
product the raw EP order statistics moments. One can also use such a result to 
compute raw moments of order statistics of a random GEP sample. 



4 Estimation and testing 

In what follows we shall discuss point and interval estimation, and also hypoth- 
esis testing inference on the parameters that index the GEP distribution. Let 
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6 = (A,/3,ay be the parameter vector. The log-likelihood £ = £(9) for a single 
observation x of the random variable X ~ GEP(A,/3, a) is 

£ = -A + logiaAfi) - alogil - e'^) + {a - \) log{l 
- e-^+'i-PW} -I3x + Ae'l^', x>0. (10) 

The score function for ([H) is given by U{e) = (d£/dA, d£/d/3, df/daV, where 

+ (g-^^ - 1)[1 - (1 - a){e'*-^'=''P(-^^"^ - 1}-'], 

e"^ - I 

X + Axe-^'[(a - i){e'i-^'=''PW - - 1], 

log(l - e-^) + log{l - e-'i+'^^'<p(-A^)}. 

As is well known, the expected value of the score function equals zero, i.e., 
E(f/(e)) = 0, which implies 

E[log{l -e-''+'^'^''P(-^^^}] 

E[(e-^^- 1){1 -(a- l)(e''-'*'=«)- 1)-!}] 
E - i)(e'-^^-p(-l^^^ _ _ i}] 

The total log-likelihood of the random sample x = (xi, . . . , XnV of size n from 
X is given by = 2"=i £^'^ and the total score function is given by U„ = 2"=i U^'\ 
where is the log-likelihood of the ith observation and U^'^ is as given above, for 
/ = 1 , . . . , n. The maximum likelihood estimator § of is obtained by numerically 
solving the nonlinear system of equations U„ = 0. It is usually more convenient, 
however, to use a nonlinear optimization algorithm (such as the quasi-Newton 
algorithm known as BFGS) to numerically maximize the log-likelihood function 

Fisher's information matrix is given by 





Ka,a 


KAfi 


I^A,o 


K = KM = n 


KAfi 








I^A,Q 







_ - 1 

dA~A 
d£_l 

da a 



log(l-e-^)--, 
a 

a 1 ^ 
- - + 1, 



(11) 



oA 



1 A 



= A-\B(X)-/3). 
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where 



1 



ae 



+ (a- l)e^B 



-4exp(-ySX) 



-j8X _ I 



1 



1 



K,^ = -j^B(X)--\-A(a-l)e''B 



X 



gA-Aexp(-/3X) _ I 



(-g^-iexp(-j8X) -1)2 



1 J a 
Ka,A = T < 1 + 



a-1 



1 -fix) _ 1 



+ (a- l)A^e^E 
1 f 1 



-2/3X -iexp(-/;X) 



- - E(X) - AE(Ze-^^) 



The above expressions depend on some expectations that can be easily solved 
computed using numerical integration. 

Under the usual regularity conditions, the asymptotic distribution of 

^fh{e-e) is N^{Q,K{er'), 

where lim„^oo = Kiff)'^. It is noteworthy that the multivariate normal 

distribution Nj,{0,Kn{6y^) can be used to construct confidence intervals for the 
parameters. In fact, an (1 - y) x 100% (0 < 7 < 1/2) asymptotic confidence 
interval for the fth parameter 61 in 6 is 

ACIi = (di - zi-y/2 V^, di + Zi-r/2 V^), 

where denotes the ith diagonal element of Kn0)~^ for / = 1, 2, 3 and Zi-y/2 is 
the 1 - y/2 standard normal quantile. 

We shall now move to hypothesis testing inference on the parameters that 
index de GEP law. Consider the partition 9 = {9j,6jy of the GEP parameter 
vector and suppose we wish to test the hypothesis "Kq '■ ^1 = ^^1''^ against the 
alternative hypothesis "Ha : 9i Of''- To that end, we can use the likelihood ratio 
(LR) test whose test statistic is given by w = 2{e(9) - £(9)}, where 9 = (9J, 9iy 
and e = ((9f^y, S^Y denote the MLEs of 9 under the null and the alternative 
hypotheses, respectively. Under the null hypothesis, w is asymptotically (as n 
00) distributed as xl^, where k is the dimension of the vector 9i of parameters of 
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interest. We reject the null hypothesis at the nominal level y (0 < y < 1) if 
w > x\ where is the 1 - 7 quantile of Using this test, one can select 
between a GEP and an EP model, which can be done by testing 'Ho : or = 1 against 
: « ^ 1. 



5 Renyi and Shannon entropies 

The entropy of a random variable X is a measure of uncertainty variation. The 
Renyi entropy is defined as /^(y) = \o%{^p {x)dx\ where y > and y ^ 1. 
We have that 

r r{x)dx = [-^^^^\ Til - e-''+'i-PW}(--i)>e-r/'-+r^»PWj;c. 

Jo ~ ^ 'Y ) JQ 

For (a - l)y + 1 > 0, we used © to expand the term {1 - e-'i+'iexp(-/3x)}(Q.-i)r 
the previous integral and obtain 

/rr w f c^'^Pe-' (-iyr((a - l)y + 1) ^ 
f{x)dx = \ — \ ) — — — :rT7^Aj{A, y), 



.00 {A(y+j)f 
■'k=0 k\(k+y) 



where Ay(^,y) = Jj^'" e-r/?A-+47+;) exp(-/3.v) ^ 

Hence, the Renyi entropy for GEP distribution can be expressed as 

= ^logLr-i y (-im(a-i)y.i) I 

Shannon's entropy is defined as E[- log /(X)] and can be obtained by taking 
limy-[-i luiy). In our case, 

E[- log/(X)] =X- \o%{aAp) + alogCl - e'^) + im{X) - XEie'l^^) 
-{a- l)E(log{l - e--^+'^-P(-A'^)}). 

After changing the variable u = e-^+^^^p(-P^) and using the expansion logx = 

-X 

T 



ZZi (valid for < X < 1), it becomes 



Efe"'®^) = — — ^— r n _ e-''+'iexp(-/3x)|Q'-l -i-2/?.v+iexp(-/Jx)j 

Jr>l po ,^ -A\k 

logKi - ^rVw + 1 = - y ^ \ 
e-^ ^ Air" Kk + a) 



(I _ g-A)a 
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Finally, using ([8]) and (fTTI) we obtain 



E[-log/(X)] = -logM/?) + logd - e-') + — - « y TTT^ 

6 Estimation of the stress-strength parameter 

We shall now discuss estimation of the stress- strength parameter for GEP dis- 
tribtuion, which is defined as 7? = Pr(y < X), where X ~ GEP{A,/3,ai) and 
Y ~ GEP{A,/3,a2) are independent random variables. Estimation of R is often- 
times of interest in statistics and has received considerable attention in the last few 
decades. Estimation of the Pr(y < X) when X and Y are normally distributed was 
considered by Govidarajulu (1967) and Church and Harris (1970). Constantine 
and Karson (1986) considered the estimation of P{Y < X), when X and Y are 
independent gamma random variables. In our case, the stress-strength parameter 
R is given by 

Q^l^P^ / (1 „-/i+iexn(-/3:c)iffi+a7-l „-/3.v+iexn(-/3.v) J.. 



(l-e-Y'^"Oo 



Let Xi, . . . ,Xn and Yi, . . . ,Ym denote independent random samples from GEP{A,/3, ori) 
and GEP(A,/3, Ofi) distributions, respectively. The total log-likelihood £r{9*) of the 
random samples Xi , . . . , X„ and Yi,. . . , Y,„ is given by 

^r{6*) = (n + m)[\og(A/3) - A] + n log ai + m log a2 - (ain + a2m) 

n 

X log(l - e-^) + (tti - 1) ^ log{l - e-''+'i'=''P(-M)} + (a2 - 1) 



(=1 

n m 



i=\ 



i=l !=1 



(12) 



The score function associated to the log-likelihood in (fT2l) is 
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Vg(er) = 0egm,des/a^,aeg/aat,dts/aa2V. where 



m 

+ ^ + - + (a2 - l)[e^-^«P(-^^.) - 1]-!}, 

1=1 

^ = - - nlog(l - e--*) + y log{l - e-i-^«P(-^-.)}, 
oa\ a\ ^ 



m 



^ = - -mlog(l y log{l -e-^^^«P(-^^')}. 

The maximum likelihood estimator of Q* is 0*, where Q* is obtained numeri- 
cally by solving the nonlinear system of equations Ur{6*) = 0. With this, we have 
that maximum likelihood estimator of is ^ = ai /(ai + a2). 



7 Applications 

In this section we fit GEP distributions to two real data sets. The first data set is 
given by Hinkley (1977) and consists of thirty sucessive values of March precipi- 
tation (in inches) in Minneapolis/St Paul. The data are 0.77, 1.74, 0.81, 1.2, 1.95, 
1.2, 0.47, 1.43, 3.37, 2.2, 3, 3.09, 1.51, 2.1, 0.52, 1.62, 1.31, 0.32, 0.59, 0.81, 
2.81, 1.87, 1.18, 1.35, 4.75, 2.48, 0.96, 1.89, 0.9, 2.05. 

The MLEs and the values of the maximized log-likelihoods for the EP and 
GEP distributions are 

A = 31.9785, j8 = 0.0186, 4p(^) = -45.7935 

and 

A = 0.8003, yS = 0.7336, a = 2.7329, £ep0) = -39.7229, 
respectively. 
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The source of the second data set is The Open University (1993). The fol- 
lowing data are the prices of the 31 different children's wooden toys on sale in a 
Suffolk craft shop in April 1991: 4.2, 1.12, 1.39, 2, 3.99, 2.15, 1.74, 5.81, 1.7, 
2.85, 0.5, 0.99, 11.5, 5.12, 0.9, 1.99, 6.24, 2.6, 3, 12.2, 7.36, 4.75, 11.59, 8.69, 
9.8, 1.85, 1.99, 1.35, 10, 0.65, 1.45. 

The parameter MLEs and the maximized log-likelihoods for the EP model are 

A = 30.8795, yS = 0.0077, 4p(0) = -75.9447. 

For the GEP model we obtain 

A = 1.9821, = 0.2369, a = 2.3144, 4p(§) = -73.6629. 

We shall now wish test the null hypothesis Hq : EP against Hi : GEP, i.e., 
Hq : a = I versus Hi : a + 1 . The values of the LR test statistic for the first and 
second data sets are (p-value: 4.9x10'"^) and 4.5636 (p-value: 3.2x10"^), 

respectively. Therefore, we reject the null hypothesis in both cases in favor of the 
GEP distribution at the usual significance levels. The plots of the fitted EP and 
GEP densities given in Figure [3] (together with the respective histograms) for the 
two data sets show that the GEP model yields better fits than the EP distribution. 




Figure 3: Fitted densities for the first (left) and second (right) data sets. 
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8 Concluding remarks 



We have generalized the exponential-Poisson (EP) distribution by defining the 
generalized exponential-Poisson (GEP) distribution. We derived important prop- 
erties of the new distribution and obtained closed-form expressions for its mo- 
ments and for order statistics moments. Applications to two real data sets were 
presented and discussed. In both applications, the GEP fit was superior to that 
obtained using the EP model. 
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